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^ Abstract 

^-H The Nordstrom- Vlasov system is a relativistic Lorentz invariant gen- 

CS| eralization of the Vlasov-Poisson system in the gravitational case. The 

CO asymptotic behavior of solutions and the non-linear stability of steady 

states are investigated. It is shown that solutions of the Nordstrom- Vlasov 
system with energy grater or equal to the mass satisfy a dispersion esti- 
f — mate in terms of the conformal energy. When the energy is smaller than 

c2 ^ the mass, we prove existence and non-linear (orbital) stability of a class of 

^ static solutions (isotropic polytropcs) against general perturbations. The 

. ^ proof of orbital stability is based on a variational problem associated to 

the minimization of the energy functional under suitable constraints. 

1 Introduction 

A classical problem in theoretical astrophysics is to establish the non-linear 
stabihty of galaxies in equilibrium. Neglecting relativistic effects and collisions 
among the stars of the galaxy, these equilbrium states can be described as 
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stationary solutions of the Vlasov-Poisson system: 

Axcfy — ^TTp^ lim 0(t,x)=O, 

|a;| — *oo 

Pit,x) = / f{t,x,p)dp. 

Here, f{t,x,p) > denotes the distribution function in phase space of the 
stars, which are assumed to have all the same mass, with t G M., x G K^, 
p € K'^ denoting time, position and momentum, respectively; (f){t, x) stands for 
the mean gravitational field generated by the stars altogether. In our units, the 
gravitational constant and the mass of each star equal one. 

A general method to approach the stability problem for an infinite dimen- 
sional dynamical system is to construct stationary solutions as minimizers of 
a suitable functional which is preserved by the evolution. The specific choice 
of the functional to minimize and of the constraints in the variational problem 
selects the type of steady states to be constructed, as well as the notion of 
distance which appears in the non-linear stability theorem. This approach was 
successfully applied to establish non-linear stability for a large class of steady 
states to the Vlasov-Poisson system, see0[IIl[T51[ni[Il[MJ[lIJIHIMl[2S]- An 
important pre-requisite for this method to yield a rigorous stability theorem is 
that sufficiently regular solutions of the dynamical system should exist for gen- 
eral initial data of the Cauchy problem (or at least for data close to the steady 
state) . In the case of Vlasov-Poisson, the existence of global classical solutions 
for general initial data has been known for some time, see [T71 [TSJ 

In this paper we study the non-linear stability for a class of stationary solu- 
tions as well as the asymptotic behaviour of general solutions to the Nordstrom- 
Vlasov system |H [S] . The latter provides a genuine relativistic generalization of 
the Vlasov-Poisson system in the following sense: It is invariant under Lorentz 
transformations and its solutions converge to solutions of Vlasov-Poisson as the 
speed of light tends to infinity. In units such that the speed of light equals one, 
the Nordstrom- Vlasov system is given by 

d^d^-A^d^^^ f f{t,x,p)^J^, (1) 



Sf - \iS^) p+il + bH-i/^V^J • Vp/ = 4/ Sq^, (2) 



where 



n— ^' S^dt+p-"^x 



are the relativistic velocity and the relativistic free transport operator. As 
compared to the Vlasov-Poisson case, the stability analysis for the Nordstrom- 
Vlasov system presents new difficulties related to the strongly nonlinear and 
hyperbolic character of the field equations. These new features make the prob- 
lem under study relevant from a mathematical point of view, besides its original 
astrophysical motivation. 
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A fundamental role in our analysis is played by the conserved quantities for 
the system ([T|-([2]). Let us mention here the energy functional, or Hamiltonian, 



■/R3 JR3 ^ JR3 ^ Jr; 



(dt^fdx 



and the Casimir functional 



Cq(/,0)= / / e'^Q{fe-^^)dpdx, 

JR3 ■/R3 



where Q : M ^ M is any sufficiently regular function. In particular, for Q{z) = 
-z'^i 9 > 1, we infer that We'^^^'^^^^'^ /Wli is constant, the case q — 1 being the 
conservation law of mass. For the purpose of the present investigation, it is 
convenient to introduce the new dynamical variable 

f{t,x,p)^e-^^fit,x,e-^p), 

in terms of which the Nordstrom- Vlasov system takes the form 



dj+ ^7p^ • V./ - V. (^V^e^^ + j . V,/ = 0, (3) 

dU-A.A^~e^^! f{t,x,p) J^^ . (4) 
The energy functional becomes 

H{f, 0, M) = £;ki„(/, + ^ / |V.0P dx+l [ {Mf dx, (5) 

^ ■/R3 ^ JR3 

where 

E^i^{l<t>)^ f I Je^^ + Wjdpdx, 

■/R3 JR3 * 

while the Casimir functional reads 



Cq(/)= / / Q{f)dpdx. 

JR3 JR3 



In particular ||/(i)||L9 is constant, for all q G [l,c»]. The main reason to adopt 
this new formulation of the Nordstrom- Vlasov system is that the field (j> does not 
appear explicitly in the definition of the Casimir functional. Let us also notice 
that the particles distribution / and the gravitational potential (f) in the energy 
functional of the Nordstrom-Vlasov system are independent variables, whereas 
for Vlasov- Poisson they are related by the Poisson formula = —p * j^- In the 

following we shall denote / simply by / to lighten the notation. 

One of the objectives of this paper is to underline the differences between 
the Vlasov-Poisson system and the Nordstrom-Vlasov system related to the 
relativistic character of the latter model. We will then focus our attention on 
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the stability analysis to a simple class of steady states, namely the isotropic 
polytropes, which in the case of the Nordstrom- Vlasov system are defined as 

/o(x,p)= (^^^^^ , E = y^e20o + |p|2. (6) 

Here fc > — 1, c > and i?o > are constants, E is the particles energy, (•) + 
denotes the positive part. The existence of a maximum Eq for the particles 
energy is necessary for the steady state to have finite energy |3]. Moreover, 
^^'0 = (I'Qi^) is the gravitational potential induced by the distribution /o and is 
a solution of the non-linear Poisson equation 

^'^" = ^'^" / / J° , „ ^P- (7) 

Our proof of stability for solutions of the form (|6l-|7| is grounded on a 
variational argument similar to the one introduced in |23pto study the orbital 
stability of polytropic spheres for the Vlasov-Poisson system and which consists 
in minimizing the energy functional subject to the constraints of given mass and 

norm of /, for some q > 1. 

This paper is organized as follows. We begin in Section |2] by proving that 
solutions of the Nordstrom- Vlasov system with energy greater or equal to the 
mass satisfy a dispersion estimate in terms of the conformal energy. This esti- 
mate does not imply that steady states solutions must have energy smaller than 
the mass, because for these solutions the conformal energy is unbounded. In 
Section [3] we state our main results on the stability of isotropic polytropes for 
the Nordstrom- Vlasov system and reduce the problem to that of minimizing the 
energy functional under suitable constraints. Some preliminary estimates neces- 
sary to solve the latter problem are given in Section [4] In particular it is shown 
that the energy of minimizers is bounded above by their mass, thus a connection 
with the result of Section |2] is established. In Section [5] we prove the existence 
of minimizers to the energy functional and show that they arise as the strong 
limit of suitable minimizing sequences. Our proof requires the minimizers to 
have energy strictly less than the mass, a property which is shown to be verified 
if the mass is sufficiently large. Finally, in Section [6] we establish uniqueness and 
spherical symmetry of the minimizer (up to a translation in space) and show 
that it is an isotropic polytropc solution of the Nordstrom- Vlasov system with 
finite radius. 

To conclude this introduction we remark that relativistic theories of gravity, 
although not physically correct, are often used as simplified models for General 
Relativity [IS]. Moreover, scalar fields play a central role in modern theories 
of classical and quantum gravity [7] . The physically correct relativistic model 
for self-gravitating coUisionless matter is the Einstein-Vlasov system, which is 
discussed for instance in [T]. Existence and finite radius property of steady 
states to the Einstein-Vlasov system have been studied in [TS] , see also [TU] , but 
the question of their stability is currently open (see however [17]). We hope 
that the present work on the stability of steady states to the Nordstrom- Vlasov 
system may contribute to a better understanding of this important problem. 
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2 A dispersion estimate 

The aim of this section is to prove a dispersion estimate for solutions of the 
Nordstrom- Vlasov system. Although this estimate will not be used in the 
following sections, it reveals an interesting link with the assumptions in our 
stability theorem, see Section |3] Define the local energy of a solution (/, 0) of 
((Sf-Q as 

the local momentum 

qi{t,x)= / Pi f dp- dt(t)di(j), 
and the local stress tensor 



, fdp+d,<i>d,<i>+h, [{Mr - (V. 



where di = d^^ ■ These quantities are related by the conservation laws 

dte + Wx ■ q = 0, dtqi + djnj = 0, (8) 

the sum over repeated indexes being understood. Upon integration, the pre- 
vious identities lead to the conservation of the total energy and of the total 
momentum: 

H{t) — edx — constant, Q{t) — qdx — constant. 

JR3 JR3 

Moreover, solutions of ([3])-(|4]) satisfy the conservation of the total rest mass: 

M — / f dpdx — constant, 
Jr3 

which is obtained by integrating the local rest mass conservation law 

■fdp. 



9tP+V,-j = 0, p = / fdp, j^e^ ( 

JR3 

We define the conformal energy as 

Sc{t) = / \x\^e{t,x)dx. 



The Nordstrom- Vlasov system is supplied with a set of initial data (/o, (/io; 
where 

/o(x,p) = /(0,x,p), 00 (a:) = 0(0, x), (/)i(a;) = 9t0(O,x). 
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In [S] it is proved that /o € , 0o € n ^ , (/ii G n launch a unique 
global classical solution of (jsJ-Q. This solution preserves energy and mass and 
has finite conformal energy for all times (provided it is bounded at time zero, 
e.g, by giving compactly supported initial data for the field). We prove the 
following 

Theorem 1 There exists a constant to > 0. depending only on bounds on the 
initial data, such that the following holds: If the initial data satisfy H > M then 

£c{t) > {H-M)t\ 

for all t > to; if H = M , then 

£c{t)>2Qot, 
again for t > to, provided Qo > 0, where 

Qo = / {x ■ q{0,x) - (j)o(f>i)dx. 



Proof: By the first of (18]) we have 



d 



-£c = 2l x-qdx, (9) 



whence, using the second equation in 

j£c = 2 [ TT{n,)dx. (10) 

■/R3 



Here Tr (rij) denotes the trace of the tensor which is given by 



It follows that one can rewrite ( [To| as 

^£c^2H + 2Q{dt<j),V^q^)-2 [ fi{t,x)dx, (11) 
dt J^3 

where /x is minus the right hand side of Q and Q is the quadratic operator 

Q{dt<P, VA) = / {{M? - (V.<^)') dx. 
By means of the identity {dtcj))'^ — dt{(j>df (j)) — 4'df(j) and using ^ we have 
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Q{dt(t),Vx(t>) ds ^ I I fj, <j) dx ds + -dt / <j) dx - j 4)o(t>idx. (12) 
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From (|9|, ([TT]) and ^ we obtain 



Scit) = £c{0) 



1 dx - 



P dx + 2Qot + Ht^ 



+2 [ [ [ ^i{(|>-l)dxdTds. (13) 

Jo Jo JR3 

Using the simple lower bound ^ — 1 > — e~'>, which holds for all ^ G M, the last 



term in ( 13 1 is bounded from below by 



Jo 



^g20 _,_ |p|2 



dpdxdrds > -Mt^ 



Substituting into (13 1 we finally obtain 



£c{t) > £c{0) - UoWh + 2Qoi +{H- M)t\ 



which yields the claim. 



□ 



Remark 1 The estimate of Theorem [T] does not apply to steady states solu- 
tions. To see this, consider /o with compact support and ^ solution of ([t]). 
The fastest decay at infinity for V(/)o is 0(|a;|~^), which is too weak to bound 
the integral \x\'^\V (f)!^]"^ dx in the conformal energy. Nevertheless the bound 
H < M will also appear as a crucial ingredient in our proof of stability for 
isotropic poly tropes. 



3 Orbital stability 

In this section we state and comment our main results on the stability of isotropic 
polytropes of the Nordstrom- Vlasov system. In order to be precise with the 
formulation of our problem, we will start by introducing some notation. For J, 
M, k positive real numbers, we denote by F^^ j the space of functions / : ^ M 
given by 

^kj = {/ e n / > a.e., ll/IUi = M, < J}. 

Moreover we denote 

D^{R^) = {0 G iioc(I^^) : V0 G and </) vanishes at infinity}, 

where the condition of vanishing at infinity means that the set {a; G : 
\4>{x)\ > a} has finite (Lebesgue) measure, for all a > 0. Functions in the space 
D-'^(]R^) satisfy the Sobolev inequality 

<?7||V0|U2, v=^7r-^^\ (14) 



7 



see [ini Thm. 8.3]. The space _D^(M^) has been extensively used in fluid me- 
chanics, see for example [TS]. 

Our first result is the existence of a minimizcr to the variational problem 

inf{i/(/,0,7/-), /eri.^j, (l)eD\ e l\ E^Uf.O)<^}, fee (0,2), 

provided the mass M is sufliciently large (depending on J and k), where H is 
the energy functional defined in Obviously, the above variational problem 
is equivalent to the following one: 

= inf{£'(/,0), / e r^^j, e d\ E^Uf.o) < 

where 

£(/, 0) = Hif, 0, 0) = £;ki„(/, 0) + ^ / I V^pdx. (15) 

By abuse of terminology, we shall continue to refer to £ as the energy functional. 

Theorem 2 For all < k < 2 and J > 0, there exists Mq £ [0,oo) such that 
l\j J < M , for all M > Mq. Moreover, for all M > AIq and any minimizing 
sequence (/„, (f>n) C F^.^ j x of the functional { 15), there exist a subsequence, 
still denoted {fn,4'n)> « sequence of spatial translations Tnfn{x,p) — /„(x + 
yn,p), Tn4>n{x) = 4>n{x + Un) , with y„ G , and a minimizcr (/o,0o) G '^a-i,j x 
such that 

\\Tnfn - /oIIli ^ 0, ||T„/„ - /o||li+i/^ -> 0, ||VT„0„ - V0o|U2 ^ 0, 
as n oo. Moreover, (/o,0o) satisfies ^ in the sense of distributions. 
Let us comment some aspects concerning the statements of Theorem [2] 

i) As in the case of Vlasov-Poisson, the use of translations in the space 
variable is necessary, otherwise starting from a minimizcr (/o,(/)o) and a 
sequence of shift vectors y„ S M.^ such that lim„^oo Ij/nl — oo, the sequence 
(Tnfa, Tn4>o) — which is still minimizing — converges weakly to zero, which 

M,J 



is not in T'^j t x . 



ii) With regard to the large mass condition, we notice first that if Mg = 0, 
then our results apply to any M > 0. If Mq > 0, then /^^ j — M, for 
all < M < Mq and in this case it is possible to construct a minimizing 
sequence of the functional £ with vanishing weak limit, sec Remark [2] in 
Section HI 

iii) In view of Theorem [l] and the previous remark, it is reasonable to conjec- 
ture that the condition I^j j < M in Theorem |2] is optimal. 

As an application of the Euler-Lagrange multipliers method we show that 
minimizers are isotropic polytropes solutions to the Nordstrom- Vlasov system. 
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Theorem 3 Let < k < 2, J > 0, M > Mq and (/o,'/'o) be a minimizer of 
15) over the space T\i j x . Then /o has the form with 



given by 



An e — , — — , C > 

" I 6 M ' M r 



MV ■ 6 7r3' y 2-A:\ Ji+i/fc r 

Moreover different minimizers differ only by a spatial translation and for any 
representative in this class the following holds: (f>o{x) ~^ 0; os \x\ —^ oo, fQ(x,p) 
is compactly supported, {fa,(t>a) is spherically symmetric with respect to some 
point in M.^ and is a time-independent mild solution of the Nordstrom- Vlasov 
system. 

Uniqueness of minimizers in Theorem[3](up to translations in space) is proved 
by showing that the Lagrange multipher Eq is the same for all minimizers. In 
the case of Vlasov-Poisson, this follows by the scaling properties of the Emden- 
Fowler equation, which is the ordinary differential equation (ODE) satisfied 
by the (spherically symmetric) gravitational potential of isotropic polytropes, 
see [13] • Here, this argument does not apply due to the strongly nonlinear 
and nonlocal character of the ODE which is the equivalent counterpart of the 
Emden-Fowler equation in the Nordstrom- Vlasov case. In fact, our strategy to 
prove uniqueness of minimizers for the Nordstrom- Vlasov system must combine 
analytical and numerical tools because some integrals involved in the associated 
ODE cannot be explicitly calculated. Our numerical/analytical computations 
reveal a monotonicity property of the set of minimizers with respect to the mass, 
which yields uniqueness by the mass constraint. 

To conclude this section, we show how the previous results yield the orbital 
stability of the minimizer solution with respect to perturbed time-dependent 
solutions of the Nordstrom- Vlasov system. We adopt the expression orbital sta- 
bility to invoke any criterium for which the orbit, described by {(T^/o, Tycj)^) ; y € 
K^}, of a stationary solution (/o, (t>o) is the set of functions which remains close 
to a perturbed Nordstrom- Vlasov solution. This concept has been widely used 
in the literature, see [53] and the references therein. 

Theorem 4 Let (/o,'/'o) be the minimizer associated to < k < 2, J > and 
M > Mq. For every e > 0, there exists 6 = 6{e) such that, for all initial data 
(/'",(/)",(/)") — (/,(/), 9((/))|(^o of the Nordstrdm-Vlasov system in the class 



and 



|ir(r,C,0i")-ff(/o,0o,o)| <5, 
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the associated solution (/, cj)) E x of satisfies, for all t > 0, 

.i^^ 11/ - '^yfoh^ + ^f, 11/ ~ ^v/oIIlw/. < e, 
inf ||V0-rj,V0o|U2 + \\dt(l)\\L^ < e. 

Proof: If the thesis of Theorem |4] were false, there would exist Eq > 0, a 
sequence (/™, '/'o"„, 0i"„) and t„ > such that 

|i^(/^<„>^«)-ff(/o>o,o)| < ^ 

andeachof ||/„(t„)-ry/o||Li, ||/n(in)-ry/o||Li+i/fc, ||V0„(i„)-rj,V0o|li2 and 
||940„(f„)||L2 is greater than Eq for all y g M'^ and n G N. Here {fm't'n) is the 
solution of ([sjl-Q associated to the initial data set (/™, (/'o"„, </>i"„)- On the other 
hand, by conservation of energy and L'^ norm, {fn{tn),4'n{tn)) is a minimizing 
sequence. Hence the above conclusion contradicts the thesis of Theorem [2j □ 



4 Bounds on the energy infimum 



Throughout the paper we denote by C any positive constant that depends only 
on M, J and k. Moreover, any subsequence of a minimizing sequence {fm4>n) is 
still denoted (/„, (/)„). We start by collecting some important estimates on the 
functions 

dp 

R 



/ fdp, A*/ = / / 



induced by an element / G V\i j. 

Lemma 1 For < k <2, let uj = 3 + k, so that l + l/u) e (6/5,4/3), and 

3 5 — w 

^-2 + 2(^731)- (16) 

Then for any f G F^j j, there exists a constant C > such that 

Pf^'^'"dx<c( [ [ f^+^l^dpdx\' (j j \p\fdpdx ' 
fi/dx<c(f f f'+'^'Y'^ ( f f \p\fdpdx 

||p/||l./5 <cf / / f^+^i'dpdx\' ([ [ \p\fdpdx^ ' 

\J«? iR3 / \JK3 Jr3 
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Proof: For all r > 1, and i? > we write 



Pf 



''dx = 



< C 



f dp+ f dp \ dx 

\p\<R J\p\>R j 



\p\fdp 



dx. 



Optimizing in R and applying Holder's inequality we obtain, for all a > 1, 



Pf 



''dx < C 



Pdp 



\p\fdp 



3q-3 
4q-3 ' 



dx 



< C 



4q-3 1 

P dp \ dx 



(4<j-3)(a-l) 

\p\fdp\ dx 



Choosing q = 1 + 1/fc, a = ui/k and r = 1 + l/w yields the first estimate on 
Pf and the one on /i/ is proved likewise. The bound on ||p/||l6/5 follows by 
interpolation with the finite mass constraint. □ 

We apply the preceeding estimates to show that the energy infimum is 
strictly positive. 

Lemma 2 For allO <k <2, J > and M > 0, I^jj > 0. 

Proof: Let (/„, (/>„) C T'Ij j x he a, minimizing sequence. As > 1 + ^, for 
all ^ G M, we have the lower bound 



Ekin{fn,4'n) > a / e^" fn dp dx + (l - a) / / \p\fndpdx 

> aM + a / 0„ / fndpdx + {I - a) I / \p\fndpdx, 



for all a e [0, 1]. Using Holder's inequality, (14) and Lemma [T] we find 

£^kin(/n,0n) > cAf - a| | 0„ 1 1 ^6 1 1 p/^J | ^6/5 + C( 1 - o) 1 1 p/„ 1 1 ^6/5 

> aM - ar/||V0„||L2||p/JL6/5 + C{1 - a)\\pfJl,„: 

where the constant C is independent from n. Letting = ||V</)„||i2 and 
Yn ~ IIp/„IIl6/5 we obtain 



> aM + -{I ~ a,^)Xl 



C(l - a) - ^arj 



Choosing a sufficiently small, precisely 



a < min 



1 C 
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we obtain £{fn, 4>n) > oAf > 0. Note that a < 1. Letting n ^ oo concludes the 
proof. □ 
The proof of Theorem 2 requires precise bounds on j. We start with a 
simple scaling argument wmch relates the energy infimum for different masses. 

Lemma 3 For all < k < 2, J > and < Mi < M2, 



r 



2 jk 



M2,J ^ Mi^^^i'-^' 

Proof: Let (/„, (/)„) C F^j^ j x be a minimizing sequence and define 
where (3 = M1/M2 and a = /3^. By direct computation, 

2~fc 



||./n||Li — M2, ll/rillii + i/'^ — 

SO that {fn,4>n) C T'li^ J X D^. Moreover 

Mo 



Mo 



ll/nl 



< J, 



whence 



4j^r< lim £(/„,0„) = ^/'= 



M2 



which concludes the proof. 



□ 



Proposition 1 The inequality I^j j < M holds, for all < k < 2, J > and 

M > 0. Moreover, the strict inequality 



'-M,.] 



< M 



(17) 



holds if the mass M is sufficiently large. 

Proof: Let q= 1 + 1 /fc > 3/2 and consider f-^ G F^^ j given by 

- J9\ ^ 



M 



X{\x\<0}{x)x{\p\<-f}{p), 



(18) 



where XA denotes the characteristic function of the set A and 



Define 60, G as 



M\ 

T 



4'aix) = -ail; 



3 

An 



1/3 



x 
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where V e C;?°(R3)^ < V < 1, V'(y) = 1, for |y| < 1, ^(y) = 0, for \y\ > 2. We 
estimate the energy of (/^, 0q) as 



dpdx+ f I \p\f-ydpdx+- f |V0apdx 



< e""M + -A/7 + a-^/JiC, 



(19) 



where, by (14), 

K = 

In particular, for a = 



\Vij\^dx > ry^^llV'llie > 
/,\j<nminf £(^,0)<M 



2/3 



5/3 



7->0 

which proves the first claim of the proposition. To prove the strict inequality 
(17) for large masses, we optimize the estimate (19) by choosing 7 = 7(0;) as 



4:K 

3M 



T 



1/6 



a. 



So doing we obtain the inequality 
where fa = f^(a) and 



(20) 



A = 



5/6 



(D 



q Iq-S 
JSg-ejVf 6,-6 . 



We choose M so large so that A ^ < 1. Precisely, 



M > 



5(9-1) 



2(9-1) 



-| l/(2q-3) 



^3(9-1) j-g 



Since A ^ < 1, we can optimize (20) by choosing a — log A. The inequality 
(20) becomes 



f(/a,0a) < A-l(l+logyl)M<M, 



which concludes the proof of (17). 



□ 



Remark 2 The constant Mq in the statement of Theorem [2] is defined as 

Mo = inf{Af' > : /^.j < M, for all M > M'}. (21) 

By Proposition [1] Mq < 00. We were not able to show that AIq = 0, i.e., that 
the strict inequality ( [l7| is always satisfied. If Mq > 0, then it follows by Lemma 
K5^that /|j J — M, for all < M < Mq- In this case, the sequence {f-y,0), with 
77 given by (18), is minimizing. Since this sequence converges weakly to zero, 
then Theorem^ is false when j = M. 
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Further information on j can be obtained by assuming that the infimum 
is achieved. For any function g, we denote by g* the non-decreasing symmetric 
rearrengement of g with respect to the x variable (see |15|). 

Lemma 4 Let {fo,(f>o) be a minimizer of the functional £ over the space j x 
D^. Then 

(i) ll/ollii+i/fc — J and (j)Q < almost everywhere; 
(a) {4>o,fa) satisfy ^ in the sense of distributions; 
(Hi) The follovjing identity is satisfied: 

\l \Wcf>o\'dx= f I .J odpdx, 

which is the relativistic analogue of the Virial Theorem; 

(iv) (/o,— ^o) '^^^^ minimizer. 

Proof: Let q = 1 + 1/fc > 3/2. Assume H/oHl^ = K < J and let K G {K, J). 
Define f{x,p) = a^fo{ax,p), 4){x) = 4)Q{ax), where a = {K / K)'^^''^''~'^\ In this 
way, ll/IUi = M, ll/IU, = k and 

£{IA)= I I Je^^o + \p\2f^dpdx 

which is a contradiction to {fo,(t>a) being a minimizer. If </)o > in a set of 
non-zero measure, then the pair {foT^\4>o\) & ^mj ^ would have energy 
strictly less than £(/o, 4>q), which is again a contradiction. For the proof of (ii), 
let C e C^(]R3). As 0o -I- e D^, for all < e M, it has to be 



: 



-/ / foCdip,x)+ / V^o-VCda;, 



which is the claim. Note that differentiation inside the integral in the kinetic 
energy is justified, since exp(j)o is bounded (by one). For the proof of the Virial 
Theorem, consider the uniparametric family of functions (/„, (pa) C F^^ j x 
given by 

fa{x,p) = /o(ax,Q:"V), (l^a{x) = (j)Q{ax) . 
Since (/o,(^o) is a minimizer, the equation 

d 



, £ {fa, (pa) 

da 



= 
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is to be satisfied, which is equivalent to identity claimed in the lemma. It remains 
to prove (iv) . It follows by the general properties of symmetric rearrengements 
that (/o,-^5) e Til J X Moreover, 



\Wcj,lfdx< / \Wcj,ofdx, (22) 
see |15l Lemma 7.17]. To reach our goal it is therefore enough to prove that 
+ |p|2/* dp dx < 



=200 



Whdpdx. 



(23) 



To this purpose we use the layer cake representation of yjl? 
we write 



For y > 



-2s 



\P\- 



ds+\p\. 



Applying this to y = —(po we obtain 
^5200 + \p\'^fo dp dx 



\/e 2s + \p\' 



fX{-<t>o<s} 



ds 



fo dp dx - 



\p\ fadpdx. 



Next we use that 



-2s 



-2s 



+ m 



-X{-4>o<s}kd{p,x) 



> 



-2s 



-2s 



-XUi<s}fod{p,x), 



see [H], eq. (3), pag. 83, and then 

I I y^e20o + \p\i fadpdx > I I 7e2^S + \p\^ fadpdx 

> I I Je-^^o + \p\^ J*dpdx. 
This concludes the proof that (/q , — '/>o) is a minimizer. 



□ 



Remark 3 Wc shall prove in Section |6] that (/o,0o) coincides with (/o,— 0o)i 
up to a translation in space. In particular, minimizers are spherically symmetric 
with respect to some point in M"^. 



5 Existence of a minimizer 

In this section we prove Theorem |2] We split the proof in several lemmas. In 
view of Lemma|4](i), it is necessary to show that the positive part of 0„ vanishes 
in the limit, which is done in the next lemma. Denote by g"*" = max(f;, 0) and 
= min((7, 0), the positive and negative part of a real valued function g, 
respectively. 
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Lemma 5 Let {fm'Pn) ^6 minimizing sequence. Then {,fn,4'n) '^^^'^ min- 
imizing sequence and 

|!V0„- V0-|U2 ^0, n^oo. (24) 

In particular, after possibly extracting a subsequence, (f>^ — > a.e. 

Proof: It is clear that {fn,4>n) C j x D^. Moreover, since |V0„p = 
\SI4>~\'^ + |V0+p, we have £{fn,(l)~) < £{fn,4>n), whence {f„,(j)~) is also a 
minimizing sequence. Now assume that ( [24| is false. Then there exists a sub- 
sequence (j)n and A > such that ||V0„ — V0~|jx,2 > A. Thus 

f(/„,0„) = i?kin(/„,0„) + ^ / (|V0-p+|V(0„-0-)|2) 

> £(/„,0;^) + 1 II V0„-V0-||i. >£(/„,(/)-) + y, 

which contradicts lim„^oo £{fn, <Pn) = hm„^oo £{fn, (t>n) = -^a/,j- 

Let (/„, 0„) be a minimizing sequence. Since (/„, (/)„) is bounded in L^+^Z'^ x 
Z?^, there exist fo S L^+^Z*^, (/ip S and a subsequence (/n,0n) such that 

/„ - /o in (^0 in i^ V</)„ - V0O in 

and so, by [121 Cor. 8.7], 

0n — *■ 0Oj pointwise almost everywhere. 

By Lemma [s] (/)o(a^) < 0, for almost all a; G M"^ and by weak convergence, 
/o > a.e. It is clear that (/o, (f>o) is our candidate for being a minimizer of the 
functional £. In the next lemma we show that the energy functional is weakly 
lower semincontinuous.. 

Lemma 6 For allO < k <2, J > and M > 0, 

lli,j>£{f\uM- 

Proof: Clearly, 

liminf / |V0„pdx > / |V(/)oprfa;. 

Moreover 

y^e^"^" + IpP \J e'^'t'o + |pp, pointwise a.e. (up to subsequences) 
and so, for all i? > and since 0o is non-positive, 

v/e20. + |p|2/„ ^ ^e^^o + \p\2f^ in ^1+1/^(5^), 
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where Br = {{x,p) € : jxp + < R'^}, It follows that 

/ / J e'^'I'o + IpI"^ f dp dx < liia'mf / We^'*" + |p|2/„ rfpdx. 
We then have 

f(/o>o) - lim / / Je^'i'o + \p\^fodpdx+l f {^M^dx 



< lim £{fnAn) = /m,j. 



which completes the proof. □ 
In the next lemma we provide a sufficient condition for the strong conver- 
gence of /„ in L^. 

Lemma 7 Let {fn,4'n) be a minimizing sequence and assume 

|1V(/.„-V0o||l^ ^0, Ti^cx). (25) 

Let M > Mq. Then 

Wfn - /oIIli ^0, as n -> c». 

Moreover, (/o, (po) is a minimizer of the functional £ over the space T\j j x 
and 

ll/n - /o||Li+i/fc ^ 0, asn^co. 

Proof: Let g = 1 + 1/fc > 3/2. To prove weak convergence in of /„, it is 
clearly enough to show that no mass is lost at infinity, i.e., that for all e > 0, 
there exists R{e) > such that 

/ fndpdx<e, 

13 J An 

where Ar = {{x,p) <eM.^ : > i? or \p\ > R}. Since 



fndpdx < / / f„dpdx+ / fndpdx 

Jm.^ J\p\>R J\x\>R.Jr^ 
< 4-^kin(/n,0n) + / / fn dp dx 

and since the term i?~-'^£'kin can be made arbitrarily small — by taking R suffi- 
ciently large — it is enough to prove that 

Ve > 0, 3R{s) > : / p/„ dpdx < e. (26) 



If (26) were false, then we could find < Q < M, a subsequence /„ and a 
sequence R{n) (depending also on Q) such that 

/ fn dp dx — Q, lim R{n) — oo. 

J\x\>R{n) JR3 
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We write 

/ /„ dp dx 

\x\>R(n) Js? 



< / / e'''"fndpdx- / / 4>nfndpdx 

l\x\>R{n) Jm.^ J\x\>R(n)JK^ 



<: / / ^e2'^" + b|2/„dpdx+ / |(/.„| / fndpdx, 

J\x\>R{n) JR^ ^ J\x\>R{n) 



where the simple bound exp(0„) — 0„ ^ 1 has been used. By (14 1 and (251, 
||0ri — 4>o\\l'^ — > 0, rt — > oo; whence, for all e > and sufficiently large n, 



[ \cj)nfdx^e. 

J\x\>R{n) 



l\x\>R(n 

Thus, by Lemma [I] and Holder's inequality. 



/ \4>n\ fndpdx ^ 

J\x\>R(n) JR3 



£(n), 



'|a;|>_R(n) 

where e(n) — > 0, as n — > oo. In conclusion 



\x\>R{n) 



^e^'f'^ + dp dx + e{n). 



We shall prove that the above inequality leads to the contradiction that I^j j 
is not the energy infimum in the space j. To this purpose, we define the 
sequence ^ ^ 

fn = af„{Px,p) X{|a;K/3-ifl(n)}, 4>n = 4>n{.Px) , 

where 

It is easy to check that the sequence (/„ , 0„) is contained in F^^ jXD^. Moreover 

£{fn,^n) = / Je^^^+\p?fndpdx+ \ j |V0„n 

\J\x\^R{n) JR3 ^ ^ JR3 J 



M-Q 
Letting n ^ oo we obtain 



\ J\x\>R(n) JM3 * y 



liminf£(/„,^„)<(^) ilt,,-Q) 
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Let us denote by F{Q) the function in the right hand side of the latter inequahty. 
It satisfies limQ^Q+ F{Q) — j, linig^^,/- F{Q) — —oo and 




for all Q G (0,M) (because of the assumption M > Mq), which leads to 
the contradiction lim„^oo •?(/«, (^n) < Im j- This concludes the proof that 
fn fo in i^- By weak convergence, |1/o||li = M and H/oHl-j < J- Thus 
{fo,<t'o) is a minimizcr by Leninia|6] By Proposition |4] ||/o||li+i/'' = J- Hence 
liminf„^oo ||/n||Li+i/fc = ll/olLi+i/': and so, after possibly extracting a subse- 
quence, /„ /, strongly in L^+^/'^ (see for instance [TS] Thm. 2.11]). Extract- 
ing a subsequence, the convergence also holds pointwise almost everywhere. A 
standard application of Egoroff's theorem shows that weak convergence in 
and almost everywhere pointwise convergence of a sequence of functions imply 
strong convergence in L^. This concludes the proof of the proposition. □ 

By virtue of Lemma [7] Theorem [2] will follow if we prove that minimizing 



sequences, after properly translated in space, satisfy (25 1. Our strategy is to 
show that it suffices to prove this for a special class of minimizing sequences, 
which enjoy some additional regularity and to which one can apply arguments 
similar to those valid for the Vlasov-Poisson system. 

Lemma 8 Let f G j being fixed such that £'kin(/: 0) < oo. 
(i) The variational problem 

has a unique minimizer ^jjf G . Moreover, ijjf is a non-positive contin- 
uous function and satisfies 

AiP^e^^ [ —J^=dp, (27) 

in the sense of distributions, 
(ii) For all (f> £ we have 

JfW-If>^-\\V(b-V^f\\h. 

Proof: Let ipn G F)^ be a minimizing sequence for the functional Jf] there 
exists tpf £ and a subsequence — still denoted ipn — such that i/;„ ^ ipf in L^, 
Vi/in — ^ Wipf in and tpn — *■ ipf pointwise almost everywhere. Moreover, by 
Lemma |5] (with kept fixed), ipj < a.e. We prove that "0/ is a minimizer by 
showing that the functional Jf is weakly lower semicontinuous. Clearly 



liminf / |Vi/;„|^da; > / \WM^dx. 
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Moreover, by convexity of the function x -^/e^^ + 

Euini^fJ) - EMn{->p7iJ) < / i'ipf-lpn) / / „ , , =fdpdx. 



R3 y/e^^f + |p|2 ■ 

The function e^'^J' /Rsi^^'''''^ + is dominated by pf and so it belongs 

to L^/^ by Lemma [ij By weak convergence, the right hand side of the last 
inequality converges to zero as n — > oo, whence 

liminf £;kin(/, V'ri) > -Ekin(/, Vy) 

n — >oo 

and the proof that ipf is a minimizer is complete. Moreover, ipf is a weak solution 



of the Euler-Lagrange equation for the functional Jj , which is ( 27 ) . By Lemma 
[l] the right hand side of (27 1 is in where j > 3/2 is given by (16 1. Whence 



ipf € Wi^^(Mr), which is continuously embedded in C(M ). In particular, ■0/ is 
bounded. Since the functional Jf is uniformly convex in n L°°, uniqueness 
of minimizers follows by standard theory of calculus of variations (see |9], for 
instance). It remains to prove (ii). For this purpose we write 

Jf{<t>) ^If = Jf{^) - JA^f) > (</> - V'/) / . , , , , „ / dpdx 

+ / Vi/'/ • (v./) - VVy) da; + ^ / \^^-\I^jfdx 
Jr^ 2 Jr3 



where the convexity of a; e^^ + |pp and the Euler-Lagrange equation ( 27 1 



have been used. This completes the proof of the lemma. □ 
Lemma 9 Let f E ^mj being fixed such that i?kin(/, 0) < oo. The equation 



(21) has a unique solution V'/ G ■ Moreover, tpf is a non-positive continuous 



function and satisfies 



e^^f(y) f f{y,p) 



\x - y\ 7r3 ^e^V-zl!/) + 



VV'f= / 1^ iLLe'^^f^y) , -'^"'"^ =dpdy, 



^ ' \x~y\ Wr3 ^e2%(-) + b|2 ^ 



f f{y,P') ^ A 
7m3 y^e^^f(y) + b'P / 
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Proof: The existence claim follows directly from Lemma [8] To prove unique- 
ness, it suffices to show that any solution of (27 1 is a minimizer of the functional 
Jf. This follows by the inequality 



Z/-JfW>^l|v^-v%||i., 



which is valid for all solutions ip G of (27) and which is derived by using 
the convexity of x — > e^^ + Now set f{x,p) = e'^'^ f {x , e"^ p) . So doing, 
(27) becomes the linear Poisson equation Aip = J^^ /(I + \p\^)~^^'^dp. Note 
that / has the same L'' regularity of /. Moreover f\p\ dp = J^s f\p\ dp < oo, 

whence, by Lemma 1 J^s /(I + e . The rest of the lemma follows 

by standard potential theory, see for instance [HI Thm. 6.21]. □ 

Corollary 1 Let (/„, </<„) C F^j j x be a minimizing sequence for the func- 
tional £ and let G he the unique weak solution of 

A^ = e^^ I —J^=dp. (28) 

Then {fm'4'fn) ^■s still a minimizing sequence and 

II W/,, - V(/)„||l2 ^0, n ^ oo. (29) 

In particular, ipf^ — ^ '/'o, as n -^ oo, pointwise almost everywhere. 

Proof: As f (/„, r/>„) > inf{Jf JV), ^ S - ^ifn^i^fj, then (/„, V'/J is still 
a minimizing sequence. Moreover, by Proposition |8] (ii) we have 

£{fn,^n)-£{fn,^fj > ^ 1 1 V0„ - V7A/„ 1 1 



and since the left hand side converges to zero, (29) is proved. Thus, after 
possibly extracting a subsequence, {ipf„ ~ 4'n) — > 0, n ^ oo, pointwise a.e. and 
since 0„ — > '/'o in the same sense, the proof of the corollary is complete. □ 
Thus the problem of proving strong convergence for V(/)„ in up to spatial 
translations has now been reduced to that of proving the same claim for the 
sequence Vy^. This problem will be addressed next, using the deep result by 
Burchard and Guo 

Lemma 10 Let {fnjipf„) be a minimizing sequence to the functional £ given by 
Corollary^ Then, there exits 4>o in such that 

i) lim„^oo ||V-i/'/„ - V0o|li2 0, 

ii) lim„_oo IIVtA/^II^s = ||V0o|li2. 
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Proof: Thanks to Lemma ji] we deduce that {fn,—'<Pf ) is also a minimizing 
sequence and verifies the inequalities ( 22 ) and ( 23 1 . By uniqueness of solutions 
to (27), see Lemma |9] is also spherically symmetric. Moreover, ipf* = 
—tjjj , by uniqueness of minimizers to the variational problem in Lemma jSj 
We now prove the first assertion of the Lemma. We notice first that since 
VV'/„ e W^^^{Br), for any ball Br = {x M.^ : \x\ < R} and the inclusion 
W^'^Br) C L'^{Br) is compact, then VV^n converges strongly in L'^{Br) (up 
to subsequences). Thus the only property that we need to prove is that 



|a:|>_R 



ivv:i 



as i? 



(30) 



where we denote Vf Integrating ( 27 1 we obtain 



J n 



from which it follows that 



\rn\r)\< 



M 



■- dpds, 



and the property (30 1 is then satisfied by the sequence ipn- conclusion, the 
sequence Vi/;* converges strongly in L^. 

The second assertion can be easily deduced from the above arguments and 
using p2l) and (p3| . □ 



Under the framework of Lemma 10 the hypotheses of Theorem 2 in [2] hold 
and as a consequence there is a sequence of spatial translations T„ such that 

lim ||r„VV'/„ - V0o|li2 =0. 

n — ^oo 

This result coincides precisely with the hypothesis of Lemma |7] and, hence, 
allows to conclude the proof of Theorem [2] 



6 Properties of minimizers 

The aim of this section is to prove Theorem [3] We start by showing that the 
minimizers constructed in the previous section are isotropic polytropes solutions 
of the Nordstrom- Vlasov system. For this purpose we use the method of the 
Lagrange multipliers. Note that we are combining some previous established 
techniques, see for example [T3J HH HH [H], together with precise change of 
scales that preserves both constraints in our minimizing problem. 

Let (/o,<^o) £ J ^ be ^ minimizer and for any e > fixed define 
S, = {{x,p) e ]R6 . g < f^(^x,p) < e-^}. Let rj e L°°{R^) be a real valued 
function with compact support such that 77 > 0, a.e. for {x,p) E M'^ysupp/o and 
supp?7C (M6\supp/o)u^e. Forte [0,T] and T = (||7?|1 1 + ||77||, + ||?7||oo)~i£/2 
we define 

ft{x,p) = a{tfM jf^l^ {a{t)x,p), 4>t{x) = (t>Q{a{t)x), 
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where 

J \\h + H\i 
M\\h + HU 

Note that /o + i^y > a.e. and, for a e small enough, 



, g=l + l/fc>3/2. 



M/2 < ll/o + triW^ <M + e/2, J/2 <\\fo + tr^Wg <J + e/2. 
From this we infer that a is a smooth function on [0,r] and 



a'{t) 



3(?-3 



a{t) 



\fo + tvh 



\\fo + tv\\: 



Moreover sup[g a"{t) is bounded. By inspection, {ft, 4>t) G T\,j j x , for all 
t e [0,r] and 



M 



ll/o 



1 EM-aUa, 4>o) 



(31) 



+ 



Mt ^ , ^ , 1 / 1 
i. . , II -Ekin ?7, (Po) + 1^ [ —7-r 
\jo + tmi 2 \a{t} 



1 



\V(f)o\'^dx. 



By a Taylor expansion at t = 0^ and straightforward estimates we obtain 

Mt 



M _ t 

a[t) 3(7 — 3 



ll/o + i^lli 



/r3 /r3 V _ /r3 /r3 /q ^ 

M Ji 



where the notation O(i^), as usual, means that the rest terms are bounded 
by Ct^, for a positive constant C depending on e, /g and ?7, but not on t. 
Substituting into (311 we get 



£{fu<i>t)-£{kAo) = t 



f f (E~Eo + cfS-')vdpdx + 0{t^), 

■/R3 JR3 ^ ^ 



where Eq and c are given in Theorem [3) Recalling the class of admissible test 
functions r; and the fact that e > is arbitrary, we conclude that E — Eq > a.e. 
on \ supp/o and /o = ( ^o-E y/(q-i) supp/o, which proves the first 

part of Theorem [3j As to the admissible range for Eq claimed in the theorem, 
we notice that, by the Virial Theorem and the definition of Eq, 



fo 
E 



dpdx 



jk 



2-k 



EqM 



The above quantity and the constant c are both positive, from which the 
admissible range of Eq is derived. 
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We proceed proving some properties of the minimizers. By Lemma [9] 

7r3 7r3 ^e^My) + |p|2 "^Ix-yl 

For all i? > 1 we split the integral in the right hand side according to — j/| < 
1/R, l/R < \x ~ y\ < R, \x ~ y\ > R. Straightforward estimates lead to 

3-7j f M 

\Mx)\<C\\iifJ,R ^ +R Pfoiy)dy+-fr, 

J\y\>\x\^R R 

which yields lim|2.|^oo 4'o{x) = 0- From this and the fact that Eq < 1, it follows 
that i/e^*^" + IpP > Eg, for \p\ > 1 or for \x\ large enough, which proves the 
compact support property of /q. Using /o = {{E — Eq)/c)\ and a change of 
variable entail 

e^-^" / , dp = Attc-'' I " Ve^ - e^'t'(Eo - E^dE < C. 

This result together with Lemma [9] implies that (po £ W^^^ , whence the char- 
acteristics of the time independent Vlasov equation are well defined curves. 
Since fo depends only on the particles energy E = •\/e^^?o~+~jpp, then it is con- 
stant along characteristics, which is the definition of mild solution. Next we 
show that (/o,0o) is spherically symmetric with respect to some point in M^. 
By using Lemma |4j (/q,— 0q) is also a minimizer. Then 



and (/o,~0o) is also a solution of Q. Writing — V'(^); ^ ~ l^^l, and 

integrating Q we obtain 

^'^r)=^ f s\y^''> ! , dpds, 

where ^' — dip/dr. From the previous equation we deduce that V^J — {x/r)^' ^ 
a.e., whence (t)Q{x) coincides with a spatial translation of (p^ix), see [3]. Since 
£^0 is the same for /o and /q , we conclude by using ^ that /o equals /g , up to 
a translation in space. 

To conclude the proof of Theorem [3] it remains to establish the uniqueness 
statement. Without loss of generality, we assume that the potential function 
is spherically symmetric with respect to the origin, i.e., 'ip{r) = (poix), r = \x\; 
recall that uniqueness is up to a spatial translation. The non- linear Poisson 
equation ^ for ip{r) is 

(rVM)' = '^'e^'^W / ^dp. (32) 

./n3 E 
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On the other hand, ^ is strictly increasing and vanishes at infinity which assures 
the existence of ro G M+ such that 



Tpiro) = log£:o- 



(33) 



Let us note that if ip{0) > logEo, then is constant in time since /o vanishes 
when e^^"^"^ > Eq. Now, we prove that tq, ^/'(ro) and il>'{rQ) are uniquely deter- 
mined by Eq for any minimizer. Let r > ro; using the Virial Theorem and the 
definition of Eq we find 



— dp ds 
E ^ 




1 

Ait 
1 

47r V 2 - fc 



— dp ds 
E ^ 



^2M-)kdpdx 
E 



EqM 



2-k 



V r > rg. 



Then, integrating (32) we have 
1 



6 



2-k 



EqM 



V r > ?'o ■ 



(34) 



Using this expression we obtain 



A7Ts^{t(;'{s)fds = 



47r 



EqM - 



k + 4: 



1 



or 



A7TS^{t(;'{s)fds = ~iP{rQ) 



E 



which allow to prove that 



-1 



47rlog£'o V2 - fc 



6 



EqM 



fc + 4 



(35) 



and 



^'(ro) =47r(log£;o)^ 



^EqM-'^ 
-fc ° 2 



M,J 



(36) 



The existence and uniqueness theory developed in [4 for ( 32 1 implies that any 
minimizer can be seen as the unique solution to ( 32 1 with initial conditions 
tp{rQ) and ip'^rg). Since these quantities are uniquely parametrized by Eq (see 
(33), ( pS] ) and ( [36| ) the question about the uniqueness of the minimizer (up to 
spatial translations) is now equivalent to prove that the spherically symmetric 
minimizer is determined by an unique Eq. 

For any Eq S (^%^J^) we have a solution ip = i/jeo to (32) with 

aere is 



initial conditions given by (33 1, ( |35| and ([36|. Associated to every ipSo ^' 
a density function / = defined by ([6]) . At this point we check that the mass 
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Compatible solutions 



Non compatible solutions 





Figure 1 : Families of solutions to ( 38 1 



associated with the functions /eo is monotone as function of Eq, which entails 
uniqueness of the minimizer by the condition ||/o||i = M. Unfortunately, these 
quantities can not be computed directly, since the initial conditions depends on 
several parameters, M, J, and k and unknown quantities like I'^j j. 



In order to skip the explicit dependence of equation ( 32 ) on Eq we scale the 
potential as 



i/'(r) = il>{hr) — \ogEa, where h ■ 



,2+fc/2 ■ 



E, 







Now ^l^{r) verifies 



E 



dp, 



with E ■ 



bP> / ^ (1 ^ E)X ^'^d initial conditions 
V^(fo) = , V''(fo) = h%p'{rf)) , fo = ro/6 . 



(37) 



(38) 



(39) 



Since the initial conditions (39 1 cannot be explicitly calculated, we consider 
numerical simulations with initial conditions at r = 0, 



^(0) = V(0) - In := a < , = 6V'(0) = 



(40) 



parametrized in terms of the new parameter a whose relation with E^ is at this 
moment unknown. Figure [T] shows these solutions for some values of a in the 
case fc = 1. Two different regimes have been detected numerically for several 
values of k. For small values of a (a < a* « 0.723 when fc = 1) the solutions 
intersect themselves; in particular, for initial conditions < a^(< a*), the 
correspoding solutions tpi and ip2 satisfy 



~1 ~2 



0' 



V''l(r~^) >V'2(^=0), 
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Initial condition and mass 




Figure 2: Masses / associated to the solutions of (38 1 



where fg, fg are the points such that V'i('^o) = V'2(?^o) = 0- However, for large 
values of a (a > 0.723 when A: = 1) the solutions mantain the order meanwhile 
they reach negative values. In particular for initial conditions < then 
ipiii^) < 4'2{f) for r < fg and 



^2 



holds. It is not hard to see that the second type of solutions is not compatible 



with conditions (39 1. In fact, by the relations (35 1, (36 1 and (39 1 one can easily 
deduce that tp'iro) must be a decreasing function of fg. The solutions with small 



values of a corresponding to conditions (39 1 arc thus the correct ones. We infer 



that the initial conditions in r = verify a strictly increasing relation with the 
values of Eq, that is. 



E'o < El e 



6 M 



^l(O) <^2(0) 



where -01 and -02 are the solutions to (38 1 with initial conditions (39 1 (deter- 
mined by £p and Eq). 

Figure [2] shows the relation of the mass in terms of the tpi- We observe that 



V'i(O) < 02(0) 



l/ill > 11/2 



Combining both estimates we conclude that the mass of the functions / is a 
decreasing function of Eq. The masses of the original / and the scaled / are 
related by 

b 



Eo 



l/ll 
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Since h/E^ is also a strictly decreasing function of Eq , we deduce that the mass 
is a decreasing function of Eq, and in consequence this proves the uniqueness of 
minimizers. 

In the remainder of this section we wish to comment on the vahdity of our 



numerical study. Let us note that the integral term in ( 38 1 can be writen 
equivalently as 



X dp - 

R3 E 



in 
k+1 



g2^(r) ^ 



k + 1 



(41) 



by using radial coordinates in the variable p and integrating by parts. In gen- 
eral these integrals cannot be computed explicitly, but this is possible for fc = 1, 
which avoids undesirable aproximation errors in the equation. However, our 



numerical simulations show that the solutions to equation ( 38 1 with initial con- 
ditions (40 1 typically have the same behavior as for fc = 1, i.e., two different 



qualitative shapes but only one of them is compatible with the constraints of 



the system (39). For fc = 1, the integral (41 1 reads 



X dp 

R3 E 



Vl-e2'AM(l - 2e2'^('')) - Se'^''''^''^ log 



l + Vl 



o2i/>(r) 



pip{r) 



In order to avoid the singularity in r = exhibited by equation ( 38 1 we have 
performed numerical simulations with initial conditions 



ip{e) = a < , ip'ie) = 0, where e = 10 



-5 



(42) 



which is a reasonable aproximation of initial conditions ( 40 1 , due to the vanish- 



ing of V'' at r = 0. The mass of the functions / can be easily computed once tp' 
is known. 

A final interesting remark about these simulations is the validation of the 
virial relation for any computed pair /, tp. Let us observe that the scaling 
deriving /, ^ from /, affects in the same way both terms in the virial relation. 
Since these functions verify 



E 



f{x,p) dx 



k+1 



f+^/\x,p)dx Vpe 



and 



\V(f>\'^dx 



dx 



R3 E 



(43) 



(44) 



where —(pix) = ip{\x\), the radial versions of these expressions allow to compute 
(by means of integrals over finite intervals) both terms. The results obtained in 
our simulations show a nearly total agreement (errors of order 10~^) between 
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both terms for each solution, indicating that these functions verify in general 
the Virial Theorem. 
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